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We study the groundstates of cold atomic gases on rotating optical lattices, as described by the Bose-Hubbard 
model in a uniform effective magnetic field. Mapping the bosons to composite fermions leads to the prediction 
of quantum Hall fluids that have no counterpart in the continuum. We construct trial wavefunctions for these 
phases, and perform numerical tests of the predictions of the composite fermion model. Our results establish the 
existence of strongly correlated phases beyond those in the continuum limit, and provide evidence for a wider 
scope of the composite fermion approach beyond its application to the lowest Landau-level. 



Ultra-cold atomic gases have become a very active field of 
study of strongly correlated quantum systems. While dilute 
Bose gases are typically in a weakly interacting regime, they 
can be driven into regimes of strong correlations by various 
means. The application of an optical lattice potential leads 
to a suppression of the kinetic energy relative to the interac- 
tion energy, and has allowed the experimental exploration of 
the quantum phase transition between Mott insulator and su- 
perfluid [1]. Rapid rotation of the atomic gas also leads to a 
quenching of the kinetic energy, into degenerate Landau lev- 
els [2], and a regime of strong interactions [3]. At low filling 
factor v (defined as the ratio of the number of particles to the 
number of vortices) this is predicted to lead to very interesting 
strongly correlated phases [4] which can be viewed as bosonic 
versions of fractional quantum Hall effect (FQHE) states [5]. 
In order to access the low filling factor regime in experiment, 
it may be favourable to exploit the strong interactions that are 
available in optical lattice systems [6, 7] for which methods 
exist in which to simulate uniform rotation (or equivalently 
a uniform magnetic field) [8-10]. This raises the interesting 
question: what are the correlated phases of atomic gases that 
are subjected both to an optical lattice and to rapid rotation? 

In this Letter, we study the interplay between the FQHE of 
bosons and the strong correlation imposed by an optical lattice 
potential. At sufficiently low particle density, the effect of the 
lattice has been shown to have negligible impact on the nature 
of the continuum Laughlin state at v = | [7, 10]. We focus 
on the possibility that there exist strongly correlated phases 
which have no counterpart in the continuum, but that appear 
as a direct consequence of both the lattice potential and rota- 
tion. To do so, we adapt the composite fermion (CF) theory 
[11, 12] which has been shown to accurately describe rotating 
atomic Bose gases in the continuum [13, 14], and apply this 
theory to rotating bosons on a lattice. Within mean-field the- 
ory, the lattice leads to the intricate Hofstadter spectrum for 
the composite fermions [15]. We predict a series of incom- 
pressible phases of bosons on the optical lattice, characterized 
by special relations of the flux density and particle den- 
sity n, and we construct trial wavefunctions describing these 
phases. From extensive exact diagonalization studies, we es- 
tablish the accuracy of the composite fermion approach, no- 
tably for states for which n = h ± ^n^; these correspond 
to incompressible quantum Hall states which have no coun- 



terpart in the continuum. To our knowledge, there has been 
no previous evidence for new FQHE states induced by a lat- 
tice potential. A previous proposal for quantum Hall states 
of bosons on the lattice [16] takes a different viewpoint, but 
remains untested. 

We study a model of bosonic atoms on a two dimensional 
square lattice and subjected to a uniform effective magnetic 
field, using the Bose-Hubbard model with Hamiltonian [8-10] 
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with aj bosonic field operators on site i, and hi = ajetj. We 
consider a uniform system with fixed average particle density 
n (per lattice site). The strength of the magnetic field is set by 
the flux density (per plaquette), defined by the condition 
that J2a Aij = 2irn^. Owing to the periodicity under — > 
+ 1, we choose < < 1. 

The single particle spectrum follows from the solution of 
Harper's equation, and takes an intricate form, often referred 
to as the Hofstadter butterfly [17]. It has a fractal structure 
consisting of q bands at rational flux density n^, = p/q. Sig- 
natures of this structure appear in the mean-field treatment of 
the Bose-Hubbard model [18, 19]. We wish to determine the 
groundstates of bosons beyond the mean-field regime, where 
interparticle repulsion leads to strongly correlated phases. We 
focus on the hard-core limit U ^> J, where the bosonic 
Hilbert-space is reduced to single occupations of lattice sites 
< rii < 1. In this limit, using standard substitutions, the 
Hamiltonian (1) can be viewed as a spin-1/2 quantum mag- 
net. The gauge fields introduce frustration, putting this in the 
class of frustrated quantum spin models where unconventional 
spin-liquid phases can appear. Indeed the Laughlin v = \ 
state studied in Ref. [10] is in the Kalmeyer-Laughlin [20] 
spin-liquid phase [31]. The strongly correlated phases that we 
describe here can be viewed as generalizations of this spin- 
liquid phase. 

Following the application of composite fermion theory for 
rotating bosons in the continuum [13], we construct composite 
fermions by attaching a single vortex to each boson. (An ex- 
plicit form for the wavefunction is presented below.) The CF 
transformation relates the flux density for the original atoms 
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FIG. 1 : Excitation gap of bosons on a lattice, with particle density n 
and flux density n$ , as predicted by a model of non-interacting com- 
posite fermions. The bright lines show parameters (n, n^) where the 
model predicts the appearance of incompressible quantum fluids, and 
include cases (where n/n^ is not constant) which are not connected 
to the continuum limit. We include data for nl = p/q, with q < 50. 

and the effective flux for CFs n*^ via 

n % = n 4> ± n > (2) 

where the two signs correspond to attaching vortices of oppo- 
site sign. Within a mean-field theory, the CFs are assumed to 
be weakly interacting, and to form a Fermi-sea which fills the 
lowest energy states of the single-particle spectrum. Incom- 
pressible states then occur when the CFs completely fill an in- 
teger number of bands. In the continuum, the single-particle 
spectrum consists of Landau-levels (LL), leading to an incom- 
pressible state each time an integer number, v* = n/n%, of CF 
Landau-levels is filled [13, 14]. Applying the same logic on 
the lattice, leads to the conclusion that the single particle spec- 
trum of the CFs is the Hofstadter butterfly [15], now at a flux 
density n*^. Owing to the fractal structure of this energy spec- 
trum, depending on n*^ there can be many such energy gaps, 
leading to many possible incompressible states. To determine 
the locations of these incompressible states, we need to know 
the particle densities n which completely fill an energy num- 
ber of bands of the spectrum of CF's at flux n*^. Generalizing 
from the continuum where the density of states in each LL is 
proportional to the flux density, an analysis of the spectrum 
leads to the conclusion that, when filling all states up to any 
given gap in the Hofstadter spectrum, the relation between n 
and remains linear [21, 22], n = v*n*^ + 5, with an offset 
5. The coefficients v* and <5 can be determined from the Hof- 
stadter spectrum by locating two points within the same gap. 
Then, using the reverse of the CF transformation (2), one ob- 
tains the lines of n, n^, on which a non-zero gap is predicted 
above the CF groundstate. 

Within a model of non-interacting CFs the relative mag- 
nitudes of gaps follow from those in the single particle CF 
spectrum. The gaps inferred under this hypothesis are shown 
in Fig. 1, in which the mode of flux attachment [correspond- 



ing to the positive and negative sign in (2)] is chosen to give 
the maximum gap. Note that the positive sign in (2) can be 
regarded either as negative flux attachment [23], or as attach- 
ment of the conjugate flux 1 — n due to the particle-hole sym- 
metry on the lattice. Indeed, Fig. 1 shows symmetries under 
ricj, «-> 1 — and n <-> 1 — n. In the hardcore limit, the 
Hamiltonian itself enjoys these symmetries, so the parameter 
space may be reduced to < n, < A. In this quadrant, 
the lines emerging from the corner with n = n$ = and 
constant filling factor v = n/n^ are the CF states expected 
in the continuum limit [13, 14]. Crucially, however, Fig. 1 
shows a large number of other lines. These correspond to new 
candidate incompressible states. 

The preceding discussion conjectures candidates for new 
types of correlated quantum liquids of bosons in optical lat- 
tices. However, given that the mean-field CF theory is an 
uncontrolled approximation, it is important to test these pre- 
dictions. There are competing condensed states on the lat- 
tice [22, 24, 25]. Even in the continuum limit, some of the 
correlated states predicted by composite fermion theory are 
replaced by other strongly correlated phases [2], with only 
v = 1/2, 2/3 and 3/4 appearing to be described in this form 
[14]. 

We have investigated the success of the CF construction for 
the Bose-Hubbard model (1) using exact diagonalisation stud- 
ies. We study the model for N particles on a square lattice 
with N s = L x x L y sites, in the presence of < < N s 
flux quanta. To limit finite size effects, we impose periodic 
boundary conditions (pbc) giving the system the geometry of 
a torus. (The nature of these boundary conditions is discussed 
further below.) 

In order to compare the exact groundstates with the CF the- 
ory it is useful to have a trial CF wavefunction. We generalize 
the continuum construction [13] to allow not only for the lat- 
tice, but also for the torus geometry, for which no convenient 
formation exists even in the continuum limit. We construct 
the trial CF state for bosons in a lattice, by writing 

«W(ri, -..,r N ) = *j({ri}) x * C F({ri}) (3) 

where "fj and ^>cf dx&fermionic wavefunctions. The factor tyj 
effects the flux attachment (2), and represents the "Jastrow" 
factor of the continuum wavefunction [13], but is now defined 
on a lattice and in the torus geometry. Instead of using the 
continuum form [26], we note that this function corresponds 
to a filled Landau-level of fermions at flux =F-/V, and we gen- 
erate 'i'j on the lattice as the Slater determinant describing TV 
fermions occupying the N lowest energy states on the lattice 
with =F_/V flux. The factor "fcF is the wavefunction of the CFs 
in the resulting effective field (2). This is described by the 
Slater determinant of the N lowest energy single-particle CF 
states $k at flux N% = ± N. For the cases derived above 
(and illustrated in Fig. 1), these numbers N and N% are such 
that the CFs fill a integer number of bands. Note that, in con- 
trast to the continuum limit where the groundstates have been 
studied within the lowest LL limit [5], our CF state (3) does 
not include a projection to the lowest LL. This is appropriate 
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FIG. 2: Average overlap of the CF states with the exact eigenstates 
in the (approximately) three-fold degenerate groundstate manifold of 
the v = | state. 

for the hard-core model that we study, since (3) vanishes when 
the positions of any two bosons coincide. 

The description of the trial state (3) is completed by dis- 
cussing the boundary conditions imposed on each of the func- 
tions. In the most general case, one introduces twisted bound- 
ary conditions for the bosons, defined by the phases 9^ = 
(9x,6y), such that magnetic translations of a boson around 
the two cycles of the torus (by in the ^-direction) act as 
\& — > exp[z^]^. In the definition of the trial state (3), one 
may choose boundary phases for the Jastrow- and CF-parts 
independently, defining 0* 9^. The sum of these phases is 
constrained to match the pbc for the bosons $ + 6^ F = 9^. 
However, this still leaves the freedom to vary 9j t — 9^. This 
freedom is a crucial ingredient to our construction: it allows 
one to generate the set of states responsible for the non-trivial 
groundstate degeneracy of these topologically ordered phases 
on a torus [27]. It is easy to show that, with this freedom, in 
the continuum limit the wavefunctions (3) reproduce the two 
continuum Laughlin states at v = 1/2 [28]. 

As an initial test, we have computed the overlaps 
(Atrial Inexact) | 2 of our trial wavefunctions (3) with the 
groundstates on the lattice at v = 1/2 as a function of n^>. 
The overlaps (not shown) are very close to those found with 
the continuum Laughlin wavefunctions [10]. Thus, there is 
little difference between the continuum [28] and lattice state 
(3), up to the flux density ~ 0.3 at which the overlap falls 
to a small value. 

Using our general construction (3), we are able to study 
for the first time the influence of the lattice structure on other 
continuum CF states. The state at v = 2/3 has a ground- 
state degeneracy c?gs = 3 [28]. We take overlaps of the CF 
trial states within the groundstate manifold composed of the 
three lowest states of the exact spectrum, and give their aver- 
age value in Fig. 2. The overlap is high and drops only above 
flux densities of ~ 0.35. Previous numerical evidence for 
the appearance of this CF state is restricted to the lowest Lan- 
dau level [5]. Our results show that, for sufficiently small n^, 
the CF state (3) also describes the groundstate for hard-core 
interactions (where LL mixing is strong). 

Let us now return to the main focus of this Letter: the new 
CF states that appear on the lattice. To investigate these states 
numerically we focus on the CF series derived from the most 
dominant gap in a subcell of the Hofstadter spectrum (cell L\ 
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TABLE I: Exact diagonalization results: gaps A and overlaps of the 
exact groundstate with the CF trial state with negative flux attach- 
ment Ocf = |{^'tiiai|*I'ex)| 2 . We also give the Hilbert-space dimen- 
sion dim(7i) for hardcore bosons. Where negative gaps are indi- 
cated, the CF state is the first excited state. Notes: ^overlap shown 
for first excited state. 



[17]), leading to a sequence with = \ — in. In order 
to be able to study several different system sizes for some 
states in this class, we select two points where (n, ru) are 
fractions with small denominators, and the particle density is 
low enough to avoid competition with the continuum Laugh- 
lin state: {\, f), and (n,n^) = (|, §). 

We find multiple pieces of evidence for the formation of 
strongly correlated incompressible phases at these values of 
(n, ritf,). First, an analysis of the eigenvalues of the single par- 
ticle density matrix of the groundstate shows that, as the sys- 
tem size N increases, there are N eigenvalues of order one. 
Thus, there is no evidence for condensation (an eigenvalue 
that grows with N), so the groundstate is likely uncondensed, 
and strongly correlated. Second, the spectra at these densi- 
ties typically show a single groundstate separated by a large 
gap (see Table I). The gaps we find are larger than the typical 
spacing of higher excited states, or the gaps at typical spec- 
tra at nearby flux densities [32]. This indicates that the sys- 
tem may be an incompressible liquid with a non-degenerate 
groundstate on the torus. This is consistent with the CF state, 
in which one expects a groundstate degeneracy of one, apply- 
ing the reasoning of Ref. [15]. 

Further direct evidence for the CF phase is obtained by tak- 
ing the overlap of the exact groundstates with the trial CF 
states (3). As detailed in Table I, we find that, in general, the 
trial CF states have significant overlap with the exact ground- 
state. Notable exceptions occur for certain cases (N = 4 for 
n = 1/7 and n = 1/9, and N = 6 for n = 1/9) where 
the exact groundstate has a different momentum from the CF 
state so the overlap vanishes identically. In these cases, we 
find large overlap of the CF state with the lowest lying excited 
state (as shown in Table I). We account for this behaviour as 
arising from the existence of a competing broken-symmetry 
"stripe" phase that is stabilized by derealization of the par- 
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tides around the short direction, similar to finite size effects 
in continuum studies on the torus [29]. This interpretation is 
confirmed by our studies at n = 1/9, which show that the 
groundstate is sensitive to the lattice geometry (two aspect ra- 
tios L x x L y = 4 x 9 and 6 x 6 are available for N = 4 at 
n = 1/9). The groundstate reverts to be of the CF form for 
the more isotropic aspect ratio. Unfortunately, no geometry 
with smaller aspect ratio is available for the systems (N = 4 
at n = 1/7 and N = 6 at n = 1/9). Still, our results indi- 
cate that, for the system at (n, n^) — (4, 1), the composite 
fermion state dominates the competing (striped) state at large 
system sizes, and maintains a very high overlap with the exact 
groundstate. A similar trend is evident for (n, n^) — , |), 
but in this case the available geometries at TV = 6 are still 
very asymmetric so we cannot confirm the preference of the 
CF state in this case. 

Overall, the values of the overlaps with the CF state are 
highly non-trivial, given the sizes of the Hilbert spaces and 
that the trial CF wavefunction has no free parameters. In 
contrast, for large system sizes, the overlap with a condensed 
(Gutzwiller) wavefunction is much smaller (< 10%), even al- 
lowing for optimization over the condensate wavefunction. It 
is clear that the CF ansatz is capturing the essential physics of 
the correlated phases. 

While a large overlap with the trial CF state is highly sug- 
gestive that the phase is of the CF type, it is very useful to have 
other tests of the qualitative features of the state. As noted 
above, the nondegenerate groundstate is consistent with the 
expected topological degeneracy of the CF state. Another im- 
portant qualitative test is provided by Chern numbers [7, 30], 
which provide a highly non-trivial test of the existence and na- 
ture of the topological order of a many-body quantum phase. 
We have evaluated the Chern number C for the groundstates 
with nonzero overlap with the CF states for systems up to 
N = 5. In all cases, we find that C = 2. This is the value 
expected for the CF phase [15]. This agreement lends very 
strong evidence that the phase appearing in the numerics is of 
the form predicted by the CF theory. 

In conclusion, we have presented numerical evidence for 
novel types of correlated quantum fluids for bosons on rotat- 
ing lattices. Our results show strong evidence that there exist 
fractional quantum Hall states beyond those in the continuum 
limit. They provide the first evidence for a wider applicability 
of the composite fermion ansatz. They also further motivate 
experimental studies of rotating gas on optical lattices, which 
would be able to probe novel aspects of the physics of quan- 
tum Hall systems. 
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